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ABSTRACT 

Context. The Rossiter-McLaughlin (RM) effect, a rotational effect in echpsing systems, provides unique insight into the relative 
orientation of stellar spin axes and orbital axes of eclipsing binary systems. 

Aims. Our aim is to develop a robust method to analyze the RM effect in an eclipsing system with two nearly equally bright 
components. This gives access to the orientation of the stellar rotation axes and may shed light on questions of binary formation and 
evolution. For example, a misalignment between the spin axes and the angular momentum of the system could bring the observed and 
theoretical apsidal motion into better agreement for some systems, including VI 143 Cyg. 

Methods. High-resolution spectra have been obtained both out of eclipse and during the primary and secondary eclipses in the 
VI 143 Cyg system, using the 0.6 m Coude Auxiliary Telescope (CAT) and the high-resolution Hamilton Echelle Spectrograph at the 
Lick Observatory. The Rossiter-McLaughlin effect is analyzed in two ways: (1) by measuring the shift of the line center of gravity 
during different phases of the eclipses and (2) by analysis of the line shape change of the rotational broadening function during 
eclipses. 

Results. We measured the projection of the stellar rotation axes using the rotation effect for both main- sequence stars in an eclipsing 
binary system. The projected axes of both stars are aligned with the orbital spin within the observational uncertainties, with the angle 
of the primary rotation axis J3p = 0.3 + 1.5°, and the angle of the secondary rotation axis /3s = -1.2 + 1.6°, thereby showing that the 
remaining difference between the theoretical and observed apsidal motion for this system is not due to a misalignment of the stellar 
rotation axes. Both methods utilized in this paper work very well, even at times when the broadening profiles of the two stars overlap. 

Key words. Stars: individual: V1143Cyg - binaries: eclipsing - Techniques: spectroscopic - Methods: data analysis 



1. Introduction 

Eclipsing binaries are great stellar laboratories for gathering in- 
formation on stellar surface structure. During eclipses, varying 
parts of the stellar disk are obscured, allowing the observer to 
gather spatially resolved information. Without eclipses, this in- 
formation is difficult to access. 

The crossing of a companion in front of a rotating star causes 
a change in the line profile of the eclipsed star as, for example, 
it first covers mainly the part of the stellar surface which is mov- 
ing towards the observer. This change in the line profile results 
in a change in the center of gravity of the line and therefore in a 
change in the measured radial-velocity of the star. The strength 
and shape of this rotation effect is a function of the projection of 
the stellar axes on the sky, its inclination (for stars with differen- 
tial rotation), the projected rotational velocity, the stellar radius, 
the radius of the companion, the stellar limb-darkening, and the 
orbital parameters of the system. 

The rotation effect was first observed by Rossiter (1924) 
in Lyrae, and by McLaughlin (1924) in the Algol system. 
The theory of the rotation effect is well understood (e.g. Kopal 
1959, Hosokawa 1953, Ohta et al. 2005, and Gimenez 2006). In 
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contrast, observations of the rotation effect in eclipsing binary 
systems are rare (e.g. Hube & Couch 1982 and Worek 1996). 
Observation and analysis of the RM effect has recently received 
renewed interest, caused by the possibility of observing the spin- 
orbital alignment for transiting exoplanet systems (e.g. Queloz 
et al. 2000 and Winn et al. 2006) and the potential to observe 
features of the planetary atmosphere (Snellen, 2004). 

For the successful observation and interpretation of the rota- 
tion effect in a planetary system, the required S/N and precision 
in radial velocity are significantly higher than those required to 
analyze the RM effect in a stellar binary system. However, the 
difficulty in analyzing a stellar binary system lies in the fact that 
one has to deal with the additional light from the eclipsing fore- 
ground star. The spectral lines of the two stars normally blend 
during the eclipses, which makes an analysis of the rotation ef- 
fect in the framework of the change of the center of gravity dur- 
ing an eclipse difficult. Nevertheless, the observation of the ro- 
tation effect is of astrophysical interest in binary systems, as it 
might reveal the orientation of the stellar rotation axes and pro- 
vide information about stellar surface velocity fields. The knowl- 
edge of these quantities might help to answer questions related 
to binary formation and evolution, and to the study of apsidal 
motion. 

Accordingly, our aim in this research is twofold: I) to de- 
velop a method for deriving information about the orientation 
of the stellar rotation axes in an eclipsing system with two 
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Fig. 1. Nonnalized spectrum of VI 143 Cyg covering the wavelength range from 5580 to 5620 A. The thin line represents the spectrum, and the 
thick line represents the convolution of the narrow-lined template with the BF obtained using the S VD algorithm. The quality of the representation 
of the spectrum by the BF and the template is also typical for other wavelength regions. 



Table 1. Parameters for VI 143 Cyg taken from ESA (1997)$, Andersen 
et al. (1987)t and Gimenez & Margrave (1985)*. Radius^ denotes the 
radius of the primary component and Radius^y the radius of the sec- 
ondary component. LJLp denotes the luminosity ratio between the sec- 
ondary and primary. 
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nearly equally bright components; II) to apply it to an astro- 
physically interesting system, VI 143 Cyg (e.g. Andersen et al. 
1987; Gimenez & Margrave 1985). VI 143 Cyg (Table 1) is a 
bright system consisting of two F5V stars, and has a high ec- 
centricity {e = 0.54) that makes it an ideal candidate for testing 
a new algorithm. Because of the high eccentricity, the spectral 
lines are not as extensively blended during eclipses. VI 143 Cyg 
is a young (2 • 10^ yr) system (Andersen et al. 1987), and the 
measured apsidal motion, i.e. the precession of the orbit in its 
own plane (0.000705+0.000041 7cycle, Gimenez & Margrave, 
1985), is only marginally compatible with what is expected the- 
oretically (0.00089±0.000157cycle, Andersen et al., 1987). The 
precession of the periastron is caused by a general relativistic 
effect and a Newtonian contribution, the latter consisting of two 
terms which are due to the deformation of the two stars by tides 
and stellar rotation. Andersen et al. (1987) suggested the pos- 



sibility that it could be possible that the tidal evolution has not 
yet achieved parallel rotation axes of the stars and the orbit. This 
would reduce the expected apsidal motion, thereby bringing it 
into a better agreement with the measured apsidal motion. 

In the following section we present our observations. 
Section 3 describes the data reduction and the two methods used 
to derive the orbital and stellar parameters. The results are dis- 
cussed in Section 4. A summary is given in Section 5. 

2. Observations 

VI 143 Cyg was observed in the summer/autumn of 2005 
and 2006 with the 0.6 m CAT telescope at the Lick obser- 
vatory, equipped with the high-resolution Hamilton Echelle 
Spectrograph. Observations of the primary eclipse 4 hours) 
were made on the night of August 29/30, 2005 with two ob- 
servations before, nine observations during, and four observa- 
tions after the eclipse. The central part of the secondary eclipse 
(^8 hours) was observed during the nights August 04/05 and 
27/28, 2005. In addition, another 20 observations were made out 
of eclipse in order to obtain an accurate orbit model. During each 
night, we also observed a set of different radial velocity standard 
stars. Before and after every exposure of VI 143 Cyg 20min) 
or a standard star, we made a Thorium- Argon (ThAr) exposure 
to obtain a wavelength scale taken close in time to the observa- 
tion, minimizing the influence of drifts in the spectrograph on 
our measurements. Data from the secondary eclipse taken on 
August 27/28 had to be discarded from the final analysis, be- 
cause of an erroneous ThAr wavelength calibration during that 
night. 



3. Analysis and results 

For our analysis we used a wavelength range from 4430 A to 
5750 A, excluding the regions with telluric lines, UJ3, and areas 
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Fig. 2. The panel shows two BFs obtained shortly before (higher radial- 
velocities), nine during, and two after the primary eclipse (lower radial- 
velocities), with the secondary BF subtracted, in velocity space. One 
can see how the part of the stellar disk which is moving towards the 
observer is initially obscured. Later in the eclipse a more central part of 
the stellar disk is covered while at the end of the eclipse parts of the disk 
that emit red- shifted light due to the rotation are covered. Note that the 
BF of the primary is moving through velocity space. These observations 
were obtained at ^ 30 min intervals. 
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Fig. 3. The same as Fig. 2 but for the secondary eclipse. One can see 
BFs of the eleven observations obtained during the secondary eclipse 
with the line profile of the primary subtracted. The change in radial- 
velocity during the secondary eclipse is smaller than during the primary 
eclipse also, the coverage of the secondary is less than the coverage of 
the primary during its eclipse. This is due to the orbital inclination not 
being 90°, and the greater distance between the two stars along the line 
of sight during secondary eclipse (See Fig. 5). 



for which an insufficient continuum correction was achieved in 
the spectra (usually the edges of the orders). The first steps of the 
data reduction consisted of bad-pixel exclusion and continuum 
normalization. The latter was carried out by fitting a polynomial 
of 5* order to the spectra and subsequently dividing the spectra 
by it. 

In the next step, a weighted mean of the two wavelength 
scales, based on the two ThAr spectra obtained directly before 
and after each observation, was constructed. The weighting took 
into account time diff'erences of the photon weighted midpoint of 
the observation and the times of the ThAr exposures. The result- 
ing error in the calculated radial-velocity due to the drift of the 
spectrograph during the observation was estimated to be 75 m/s 
on the basis of comparisons of several ThAr exposure pairs. 

The Broadening Function (BF) was calculated to retrieve the 
line profiles of the two stars, governed by the velocity fields on 
the stellar surface and stellar rotation. The BF represents the 
function that projects a narrow-lined template onto the broad- 
ened spectrum. We calculated the BF using the singular value 
decomposition (SVD) technique. This approach is clearly de- 
scribed in Rucinski (1999), whereas the SVD algorithm itself 
can be found in Press et al. (1992).^ In order to calculate the 
BF one needs a narrow-lined template; we used the spectrum of 
the F7V star HD 222368 (Udry et al. 1999), obtained as one of 
the radial- velocity standard stars. Before using the spectrum of 
HD 222368 as a narrow-lined template we deconvolved it using 



^ We took all singular values until the first derivation of the difference 
between the observed spectrum and the template broadened by the BF 
with respect to the singular values used for the calculation of the BF 
effectively reached zero. 



the maximum-likelihood method in conjunction with a kernel 
calculated with the same code as applied in Section 3.2. An iter- 
ative approach was used to identify the kernel that gave the best 
results. The most satisfactory results are achieved when using 
a kernel with a projected rotational velocity (ysini) of 5 km/s. 
Fig. 1 shows a part of a typical spectrum obtained from the 
VI 143 Cyg system. It also shows for comparison the convolu- 
tion of the narrow-lined template with the BF. As VI 143 Cyg is 
a double-lined binary system with two components of the same 
spectral type, the BF itself consists of two peaks that represent 
the broadening functions of the two stars, shaped by the corre- 
sponding rotation and velocity fields. 

In this study we are primarily interested in the data taken 
during the eclipses and their interpretation in the framework of 
the RM eff'ect. Because of the high eccentricity {e = 0.54) of 
the system, the midpoints of the eclipses do not occur when the 
radial-velocities of the stars are equal, but shortly before and 
after the time of equal velocity, for the primary and secondary 
eclipses, respectively. During the middle and the end of the pri- 
mary eclipse the lines of the two stars are blended. The same 
situation occurs during the beginning and the middle of the sec- 
ondary eclipse. Hence, although the line blending is not as severe 
as in the case of a circular orbit, the light from the foreground 
star cannot be ignored. 

We followed two diff'erent approaches to derive the spin axis 
of the eclipsed star with blended lines. The first method is ex- 
plained in Section 3.1. We first show how the influence of the 
foreground star is subtracted, and subsequently how the cen- 
ter of gravity from the BF of the eclipsed star is calculated. In 
Section 3.2 the BFs of both stars are used. Here we do not use 
the center of gravity of the measured BFs to derive the RM ef- 
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Fig. 4. The radial- velocity measurements of the primary and secondary components of 
VI 143 Cyg, and the orbital solution, are plotted against orbital phase. The lower panel 
shows the difference between the best fit and the actual measurements of the radial- 
velocities of the primary and secondary component. The midpoint of the primary eclipse 
occurs at a phase of ^ 0.03, and the midpoint of the secondary eclipse at a phase of ^ 
0.77. Please note that the systemic velocity in this graph and in Figs. 6 and 7 is already 
subtracted 



Fig. 5. The orbit of the binary system 
VI 143 Cyg shown from above the orbital 
plane. The solid line represents the orbit of 
the primary component and the dashed line 
the orbit of the secondary. The lines from the 
center of gravity towards the orbits indicate the 
position of the periastron. The big dots indicate 
the positions of the stars at time of mid primary 
eclipse. 



feet, but the shape of the BFs and their change during the eclipse. 
The measured BFs are compared to simulated BFs, and thereby 
the parameters that govern the rotation effect are derived. 

3.1. Method 1: The BF's center 

Before determining the parameters involved in the RM effect, an 
orbital model has to be obtained. To derive the radial- velocities 
of the components out of eclipse, we fitted two Gaussians to 
the two peaks in the BFs. A fit was applied to extract the 
orbital parameters. In this work, we adopted the orbital pe- 
riod (7?6407568 ± 6 x 10"^) for all fits, since it is derived at 
much higher accuracy from eclipse photometry than from radial- 
velocity variations (Gimenez & Margrave 1985). The inclina- 
tion of the orbit (87.0 ±0.1°), and the sizes of the components 
1.346 ± 0.0237?o and 1.323 ± 0.0237?© have been adopted from 
Andersen et al. (1987); these are needed for later analysis. The 
fitted parameters are shown in the second column of Table 2 with 
their l-cr uncertainties. In addition, the orbital parameters given 
by Andersen et al. (1987) and Gimenez & Margrave (1985) are 
shown for comparison in column five. 

The tomography algorithm of Bagnuolo & Gies (1991) is 
used to disentangle the primary and secondary spectra. This al- 
gorithm uses, the spectra obtained at different phases of the or- 
bit, and the orbital parameters of the system, as input. It starts 
with two synthetic spectra without spectral lines for the two 
components in VI 143 Cyg. For all observations taken outside 
the eclipses, it shifts the observed spectra in the rest- frames of 
each component using the newly obtained orbital parameters. 
Subsequently, the synthetic spectra are compared with the ob- 



served spectra. The mean of all the differences between the syn- 
thetic and observed spectra is added to the synthetic spectra. This 
complete process is repeated 50 times, but in our case it con- 
verged after only a few iterations. 

In the next step, the spectrum of the foreground star is sub- 
tracted. For an observation out of eclipse this is straightforward, 
using the spectrum shifted in velocity space to the appropriate 
position. During eclipses, one has to incorporate the change in 
the light ratio of the two stars due to the eclipses. For this we 
assumed a linear limb-darkening law with a limb-darkening co- 
efficient (u) of 0.6 for both stars. Subsequently, the BF was cal- 
culated with only one star in the spectrum. Figs. 2 and 3 show 
the BFs of the primary and secondary stars during their eclipses, 
after subtraction of the foreground star. Note that, for computa- 
tional reasons, the continua of the observed spectra were set to 
zero and the signs of the spectra have been changed; this results 
in positive BFs. 

The center of gravity of the absorption lines can now be ex- 
tracted from the BFs and used to calculate the radial-velocity 
of the eclipsed star, including the radial- velocity anomaly intro- 
duced by the RM effect. Using the orbital model and the formula 
by Kopal (1959) and Hosokawa (1953)^ to calculate the rotation 
anomaly, the orbital parameters and the parameters relevant for 
the RM effect of the two stars can be derived jointly with a x'^ 
minimization algorithm; therefore vsini and the angles jSp and fis 
enter as new parameters. JS is the angle between the stellar spin 
axis projected onto the plane of the sky and the orbital spin axis 
projected onto the plane of the sky. 13 = 0° would indicate a pro- 
jected rotation axis perpendicular to the orbital plane, whereas 



^ Equation 5 has to be adopted for eccentric orbits. 
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Fig. 6. The rotation effect during the primary eclipse. The shift of the 
center of gravity is plotted against the orbital phase. The upper panel 
shows the measured radial velocities along with the best fit. In the sec- 
ond panel the radial velocity due to the orbital motion is subtracted from 
the data and the fit. The residuals between data and fit are shown in the 
third panel. 



J3 = 90° would indicate that the rotation axis lies in the orbital 
plane. The longitude of the ascending node of the orbit (Q) is 
not known, resulting in an ambiguity in the sign of the angle jS. 
In our definition, a positive JS indicates that the RM effect, inte- 
grated over the complete eclipse, would give a positive residual 
in radial-velocity. In that case, the companion spends a longer 
time in front of the stellar surface moving towards us, than that 
part of the stellar surface that is moving away from us. The lin- 
ear limb-darkening coefficients of both stars are fixed during the 
fits to 0.6 as they are only weakly constrained by our fits. In the 
fit, all data points out of eclipse and all data points during the 
primary and secondary eclipses are fitted simultaneously. Our 
radial- velocity data points and the best fit can be seen in Fig. 4 
for the complete orbit including the two eclipses. Fig. 5 shows 
the orbit from a point above the orbital plane. Figs. 6 and 7 dis- 
play the radial- velocities of the primary and secondary stars dur- 
ing and around their eclipses. The best fit parameters are listed 
in the third column of Table 2 with their 1-cr uncertainties, as 
derived from the fit- AH radial-velocity measurements are 
given with their epochs in Table A. 1 in the Appendix. The radial- 
velocities directly after the primary eclipse, and the data point at 
a phase of 0.7, have been extracted using a Gaussian fit to the BF 
of one star after the subtraction of the other star in the spectrum; 
however, there is no significant change in the derived parame- 
ters if these data points are omitted. The average uncertainties 
in the radial- velocities, out of the eclipses, are 0.15 km/s for the 
primary and 0.25 km/s for the secondary. The uncertainties of 
those radial-velocity measurements that were determined from 
the centers of the lines, we estimated to be 0.20 km/s for the pri- 
mary and 0.32 km/s for the secondary. This was calculated by 
comparing them with the velocities obtained from a Gaussian fit 
outside the eclipses. 

3.2. Method 2: Variation of the BF profile 

In our second approach, to derive information about the binary 
orbit and the orientation of the rotation axes, we simulated the 
shape of the BF of the two stars, as governed by the orbital mo- 
tion, stellar rotation, orientation of the stellar spin axes, velocity 



-20 






-25 
-30 






: 

r 




-35 






5^-40 


: 


^^^^ ^ : 


5-45 


- 


- 








•"S 2 
o 

^ 


: 








— _ 9 
P 


_ 


^^^"^ - 


o ^ 






a: 1 .0 






0.5 
0.0 
-0.5 




^ ^ -- -- ^ ■■ ^ ^ 






-1.0 







0.75 0.76 0.77 0.78 0.79 

Phase 

Fig. 7. The same as for Fig. 6 but this time for the secondary eclipse. 

fields on the stellar surface, limb-darkening and, in case of a 
measurement taken during an eclipse, the fraction of the stellar 
disk which is covered by the companion at the time of the mea- 
surement. These simulated BFs were subsequently compared to 
the measured BFs, and the relevant parameters determined. With 
this method, we not only use the first moment of the stellar ab- 
sorption lines during the analysis, but we utilize the complete BF 
and its change over the course of the eclipse as a diagnostic tool 
for determining the orbital and stellar parameters. 

We simulated the rotation profile of the two stars in the 
VI 143 Cyg system with a few thousand elements across the vis- 
ible half-spheres with equal surface brightness. We included lin- 
ear limb-darkening, asynmietric macro-turbulence and solar-like 
diflferential rotation (See Gray 2005). Accordingly, besides the 
parameters which were fitted in Section 3.1, another parameter, 
the Gaussian width of the macro-turbulence, is required. We kept 
the cr of the Gaussian for the tangential and radial- velocity fields 
^RT and their covered surface fraction equal, as the quality of the 
fit did not improve by including them as free parameters. The in- 
clination of the orbit and the sizes of the two stars can be varied 
in our fits. However, in the final fits presented here, we used a 
linear limb-darkening coefficient of 0.6, solid body rotation, and 
the inclination and radii from the literature (see also Section 4). 

Each exposure of VI 143 Cyg had a duration of ^20min. 
Provided that one takes the photon midpoint of the observation, 
this poses no serious problem for the analysis of the center of 
gravity of the lines. However, looking at the line profile, one 
has to take two eflTects into account. First, the radial- velocities 
of the two stars change during a 20 min exposure. This eflfect is 
strongest in our data set at the time of the primary eclipse, where 
the change in radial- velocity is » 3 km/s during one exposure. 
This artificially widens a BF taken during a 20 min exposure 
relative to a BF which would have been taken instantaneously. 
During the eclipses, a second eflTect becomes important; the cov- 
erage of the eclipsed stars can change considerably during 20- 
min. Again, this eflfect is stronger during the primary eclipse than 
during the secondary eclipse, as the primary eclipse is deeper 
and shorter. Sine we need an exposure over several minutes to 
receive a sufficiently high S/N ratio^, we introduced the same 
"smearing" in the simulated BFs. We calculated three BFs for 
each exposure and stacked them together: one simulated obser- 
vation 6.6 minutes before the photon midpoint, one at the pho- 

^ Out of eclipse the spectra have a S/N between 40 and 50. 
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Fig. 8. BFs during the primary eclipse centered in the frame of the primary. Each panel represents one measurement. The solid line represents 
the measured BF, the dashed line shows the best fit BP, while the dotted line represents the diff'erence between the fit and the measurement. 
The numbers in the upper right corner of each observation indicate the time diff'erence between the photon midpoint of the observation and the 
midpoint of the eclipse in hours. In the upper left corner one can see the uncovered part of the primary half-sphere for the photon midpoint of 
that measurement if the star was ecUpsed during this measurement. One can see that the BF of the secondary moves closer to the primary BF in 
velocity space over the course of the eclipse. Please note that the scale of the ordinate is arbitrary. The panel in the lower right corner shows the 
BF of the primary for an observation out of eclipse. 



ton midpoint and one 6.6 minutes after the photon midpoint. We 
convolved each simulated BF with a Gaussian of cr = 2.5 km/s, 
representing roughly the resolution of the spectrograph in the 
wavelength range used. 



In addition to the parameters describing the binary system, 
one parameter is included that scales all simulated BFs to the 
observed BFs. The primary star has a slightly higher luminosity 
than the secondary. Also, the template used might fit one of the 
two stars better than the other one. Therefore, we included an- 
other parameter scaling the height of the kernels of the two stars 
relative to each other. As mentioned in the beginning of this sec- 
tion, all spectra have been normalized. During the eclipses the 
absolute amount of light changes. Hence the depths of the ab- 
sorption lines change relative to the normalized continuum, and 
therefore the height of the BFs also change. This has to be incor- 
porated in the calculations of the BFs. We performed a fit using 
all 46 spectra obtained out of eclipse and during the primary and 
secondary eclipses. The derived values for the parameters are 
given in column four of Table 2. The measured BFs and the best 
fits can be seen in Fig. 8 for the primary eclipse and in Fig. 9 for 
the secondary eclipse. The uncertainties were calculated using 
the bootstrap method described in Press et al. (1992). 



4. Discussion 

4.1. Orbital parameters 

Using the methods described in Section 3.1 and Section 3.2 we 
obtained values for the orbital parameters of VI 143 Cyg and for 
the stellar parameters of the two stars. The time of periastron 
passage of the primary (T), the longitude of the periastron (oj), 
and the eccentricity (e) were determined in all three fits, and 
agree with each other to within the l-cr level. Using the val- 
ues from Gimenez & Margrave (1985) for oj (48.26 ± 0.01°) 
and the apsidal motion rate (0.000705+0.000041 Vcycle), and 
taking the time diflTerence into account, one would derive an 
0) of 49.31° ± 0.06° for the time of our observations. This is 
also the value stated in Table 2. The parameter for the semi- 
amplitude of the secondary (K^) (90.0 ±0.1 km/s) falls outside 
the l-cr range of the literature value (91.1 ± 0.4 km/s). This dif- 
ference can also be seen in the values for the projected semi- 
major axis of the system (a sin /). The uncertainty in the radial- 
velocity of the center of gravity of the system (y) includes the un- 
certainty of the radial- velocity in our template star, HD222368 
(5.6 ± 0.3 km/s, Udry et al. 1999). Our calculated masses for 
the two components in VI 143 Cyg are Mp = 1.355+0.004 M© 
and Ms = 1.327+0.003 M©. These values lie in between the val- 
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Fig. 9. The same as Fig. 8 but for the secondary ecUpse. One can see that the BF of the primary moves further away from the secondary BF in 
velocity space over the course of the eclipse. The panel in the lower right comer shows the BF of the secondary for an observation out of eclipse. 



ues calculated by Andersen et al. (1987) = 1.391+0.016 M© 
and Ms = 1.347+0.013 M©, and the values given by Gimenez & 
Margrave (1985) Mp = 1.33+0.03 and = 1.29+0.03 Mq. 



4.2. Stellar parameters 

The stellar parameters obtained in Section 3.1 are derived by 
analyzing the shape of the rotation anomaly in the radial veloc- 
ity. The method relies on a clean subtraction of the foreground 
star from the obtained spectra before calculating the BF's, which 
depend on orbital parameters derived from out-of-eclipse mea- 
surements, values for the stellar radii, the stellar limb-darkening 
and the light ratio between the two stars during the subtraction 
process, which have been taken from the literature. 

In Section 3.2, the change in the shape of the absorption lines 
is used instead of the rotation anomaly. Looking at Figs. 8 and 
9, one can see that the simulated BFs are somewhat "rounder" 
during eclipses than the observed BFs. This can clearly be seen 
during the central phase of the primary eclipse (Fig. 8 panels in 
the second row) and during the secondary eclipse (Fig. 9). 

The agreement between the measured BF and the observed 
BF can be improved if u is included in the fit. The derived val- 
ues for u would be 0.9 + 0.1 for the primary and 0.8 + 0.1 
for the secondary. As these values are probably too high (Gray 
2005) and the derived values for /3 would not change signifi- 
cantly (fip = -0.6 + 1.4° and fis = -0.3 + 1.2°), we kept u fixed 
to 0.6 in the final fits. 



Including solar-like diff'erential rotation, the orbital inclina- 
tion and the stellar radii as free parameters in the fits also leads 
to a better agreement between data and simulation. However, 
the derived diff'erential rotation parameters are negative for both 
stars; the angular rotation speed is faster at the poles than at the 
equator. Furthermore, the fitted radii are not in agreement with 
the literature values; the primary radius is increased relative to 
the literature value and the secondary decreased relative to its lit- 
erature value. The value for the orbital inclination of VI 143 Cyg 
in the literature is reproduced by our fit, but weakly constrained. 
A negative diff'erential rotation parameter would give the BFs 
during the eclipses a more box-like shape. A bigger diff'erence 
between the stellar radii would further improve the agreement 
between simulation and data during the central part of the pri- 
mary eclipse. However, a bigger diflference in stellar radii would 
also mean a bigger diff'erence in stellar masses, which would not 
agree with the nearly equal semi-amplitudes found in the radial- 
velocity. 

We favor a different explanation for the diff'erence of "round- 
ness" between the measured and simulated BFs. Looking at an 
absorption line originating from one small area of the stellar disk 
(i.e. no rotational broadening), the core of an absorption line is 
formed higher up in the stellar atmosphere than the wings (see 
Gray 2005), as the same optical depth is reached earlier in the 
core than in the wings of the absorption line. One can now com- 
pare the eflfect that limb-darkening has on the core and the wings 
of the absorption lines. In the core of an absorption line, a high 
optical depth is reached in a high layer of the atmosphere look- 
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Table 2. Derived parameters of V1143 Cyg are given together with their formal errors. The reduced of the orbit fit is 0.95 and the reduced 
of the joint fit is 0.96. The average uncertainty in radial-velocity for the primary is 0.16km/s and 0.27 km/s for the secondary, in the joint fit. For 
comparison, values given by Andersen et al. (1987) and Gimenez & Margrave (1985) are also shown. The second and third columns present the 
parameters obtained in Section 3.1. The fourth column shows the parameters derived with the method presented in Section 3.2. For the reasons 
mentioned in Section 4 we consider the values found in Section 3.2 to be our best parameters and the formal errors given for v sin / and ^rt as too 
small. 



Parameter 


Center 


Shape 


Andersen et al. (1987)t 




Orbit 


Joint fit 


Gimenez & Margrave (1985)* 


T [JD-2400000] 


53536.130+0.002 


53536.131+0.002 


53536.1317+0.0006 




Kp [km/s] 


88.1+0.04 


88.1+0.1 


88.01+0.05 


88.2+0.2t 


K rkm/<jl 
r^s LKiu/hj 


90.1+0.08 


90.1+0.2 


89.9+0.1 


91.1+0.4t 


e 


0.538+0.001 


0.538+0.001 


0.5378+0.0003 


0.540+0.003t 


con 


49.1+0.2 


49.1+0.2 


49.27+0.05 


49.31+0.06* 


a sin / [Rq] 


22.67±0.03 


22.67+0.03 


22.64+0.02 


22.78±0.08t 


7 [km/s] 


-16.8+0.3 


-16.8+0.3 


-16.8+0.3 


-16.5+0.7t 


Mp sin^ / [Mo] 


1.357+0.005 


1.357+0.008 


1.350+0.004 


1.386+0.016t 


M, sin^ / [Mq] 


1.327+0.004 


1.327+0.007 


1.322+0.003 


1.341+0.013t 


vsin/p [km/s] 




16.9+1.0 


19.6+0.1 


18+2t 


vsin/^ [km/s] 




28.0+5.0 


28.2+0.1 


27+3t 








3.4+0.1 




(rtS [km/s] 






3.3+0.1 








0.5+4.0 


0.3+1.5 








-3.9+4.0 


-1.2+1.6 





ing down on a stellar atmosphere. This also means that looking 
from an angle into the stellar atmosphere (i.e. at the limb) will 
result in a small change in the average height at which the core of 
the line is formed: Only the outer cooler parts of the atmosphere 
are probed. For the wings of the absorption lines the situation is 
diflTerent; the height at which the wings of the observed absorp- 
tion line are formed changes more from the center to the limb of 
the stellar disk. Therefore, in going from the center to the limb 
of a stellar disk the absorption lines become more box-like. The 
particular shape change of an absorption line originating under 
diflferent angles relative to the observer depends on the details of 
the stellar atmosphere. This would result in a non-linear limb- 
darkening law. Out of eclipse, the light is integrated over the 
whole stellar surface, but for example, during the central part of 
the primary eclipse light is only received from the outer part of 
the stellar disk (See Fig. 8, icon in the panel of the third column 
and second row). Therefore, the shape change of the absorption 
lines originating from diflTerent parts of the stellar disk might ex- 
plain the diflTerence between the observed BFs and the simulated 
BFs (e.g., Pierce & Slaughter 1982, Balthasar 1988 or Hadrava 
2007). 

Using the broadening of about 1000 lines, we obtained 
macro-turbulence parameters (^rt) of 3.4 km/s and 3.3 km/s 
for the primary and secondary, respectively. For comparison, 
the macro-turbulence parameter found for the Sun, as a disk- 
integrated star, for weak and moderately strong lines is ^ 
4.0 km/s (Takeda, 1995). 

The values obtained for the vsini of the secondary, using the 
two different methods described above, agree with each other 
and are also consistent with the literature value. Note that there 
might be a systematic diflTerence between the obtained radial- 
velocities of the secondary during the secondary eclipse and the 
fit to the radial-velocities (Fig. 7). This might be caused by too 
high an uncertainty in the orbital parameters used during the data 
analysis in Section 3.1. The values obtained for the primary vsini 
do agree with the literature value within their 1-cr error, but do 
not agree with each other within their l-cr errors. The method 



applied in Section 3.1 uses the amplitude of the RM eflTect to de- 
rive the value for the projected rotational velocity. It assumes a 
linear limb-darkening coeflficient. As mentioned in the last para- 
graph, this might not be sufl&cient. Also in Section 3.2 a linear 
limb-darkening coefl&cient is used, but vsini is derived not only 
by the use of the RIM eflTect, but also by the shape of the BF out- 
side of the eclipse. Therefore, we disregard the value for vsini 
obtained in Section 3.1 and take the value of 19.6 km/s from 
Section 3.2 as the projected rotational velocity for the primary 
star in VI 143 Cyg. Due to the remaining mismatch between the 
data and the simulation, which might be due to the use of too 
simple a model, we consider the formal errors for ^rt and vsini 
as too small. 

We would like to point out that the methods used are still in 
their infancy. Even with the derived values for the stellar radii, 
the orbital inclination, limb-darkening and diflTerential rotation 
are not free of errors yet; these parameters are not normally ac- 
cessible, or only with great diflficulty, via spectroscopic data. 



4.3. Orientation of the rotation axes 

The main focus of this work is the robust determination of the 
orientation of the stellar rotation axes. The values derived for the 
projection of the rotation axes onto the plane of the sky are given 
in the last two rows of Table 2. The J3s derived from the two dif- 
ferent methods agree to within their errors despite the diflTerent 
systematic problems of each method. Therefore, the projections 
of the rotation axes of both stars are, to within their uncertain- 
ties, perpendicular to a vector that lies in the plane of the orbit 
and is perpendicular to the line of sight. What does that mean 
for the orientation axes of the stars? It is highly unlikely that the 
geometry is such that we see the projection of the rotation axes 
perpendicular to the orbital plane and the stellar rotation axes 
have a high inclination towards the observer. We therefore con- 
clude that the stellar rotation axes are normal to the orbital plane, 
and aligned with each other and the rotation axis of the system. 
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This leaves the theoretical prediction of the apsidal motion 
unaltered for VI 143 Cyg. Hence, the difference between ex- 
pected (0.00089±0.00015°/cycle) and measured apsidal motion 
(0.000705±0.0000417cycle), which just Hes outside the l-cr er- 
ror bars, is also unchanged. The effect of a misalignment be- 
tween the stellar rotation axes and the orbital spin axis on the ap- 
sidal motion has been studied by Kopal (1978), Shakura (1985), 
Company et al. (1988), and Petrova & Orlov (2003). The con- 
tribution of the stellar rotation to the advance of the longitude 
of the periastron is reduced if the stellar rotation axis is tilted 
against the orbit spin axis until finally, when the axis of stellar 
rotation lays in the orbital plane, its contribution is half as large 
and with the opposite sign as when the stellar and orbital axes 
would be parallel. In this situation, it contributes to a retrograde 
rotation of the periastron. The contribution of the stellar rota- 
tion to the apsidal motion depends not only on the orientation 
of the axis, but also on the square of the angular stellar rotation 
rate. As one measures only v sin /, a greater inclination towards 
the observer would mean a higher angular stellar rotation rate, 
and therefore a greater contribution of the rotation term to the 
overall apsidal motion. We calculate that if the rotation axes of 
both stars would lie in the orbital plane, but have no inclination 
towards the observer, then the complete apsidal motion would 
be 0.00073°/cycle. If the rotation axes would have an inclination 
towards the observer of / ^ 70° in either of the two stars, or of 
/ » 60° in both stars, only then would the expected and measured 
apsidal motion be in agreement. The secondary, due to its higher 
V sin /, has a larger influence on the rotational term of the apsidal 
motion than the primary component. As already pointed out, it 
is very unlikely that the stellar and orbital spin axes span a large 
angle, while their projections on the sky are, in their uncertain- 
ties, parallel. 

Our findings do not support the hypothesis advocated by 
Petrova & Orlov (2003), that a misalignment of the stellar ro- 
tation axes with the orbital spin could bring the theoretical and 
measured apsidal motion for a number of binary systems, includ- 
ing VI 143 Cyg, into better agreement. 

Our work has excluded the option of a misalignment between 
the stellar rotation axes as a possible explanation for the dif- 
ference between the expected and measured apsidal motion in 
VI 143 Cyg. It is therefore interesting to look at other possibil- 
ities that might explain this diff'erence. As the apsidal motion 
constant (^2) is an important source of uncertainty in the calcu- 
lation of the expected apsidal motion, a new calculation of the 
apsidal motion constant for VI 143 Cyg using modern codes for 
stellar evolution might be of value. 

In our analysis of the orbital data we found no indication of a 
third body in the VI 143 Cyg system, whose influence might also 
alter the apsidal motion. However, because of the limited cover- 
age in time (one year) and the limited accuracy in radial-velocity, 
combined with the possibility that the orbit of a third body could 
have a lower inclination, a third body cannot be excluded. 

The alignment of the stellar rotation axes could also set a 
lower limit to the age of the system, if the axes were not aligned 
at the birth of the system. However, this is not a straightforward 
argument, since, during the pre-main- sequence phase, synchro- 
nizing forces must have been larger, due to the larger sizes of 
the stars. Proper modeling of the evolution of VI 143 Cyg might 
reveal whether it had enough time to align its axes. If the time- 
span to align the axes is longer than the lifetime of the system, it 
would mean that VI 143 Cyg, with its high eccentricity of 0.54, 
was born in this way, with all spin and orbital axes aligned. 



5. Conclusions 

We measured the Rossiter-McLaughlin efl'ect for both compo- 
nents in the binary system VI 143 Cyg. We developed two diff'er- 
ent methods to derive the angle /3 between the stellar spin axis 
projected onto the plane of the sky and the orbital spin axis pro- 
jected onto the plane of the sky. 

Using the first method, i.e. by determination of the center of 
the broadening function, we showed how it is possible to subtract 
the absorption lines of the foreground star during an eclipse from 
the spectrum. This made it possible to use the shift of the cen- 
ter of gravity of the absorption lines as a proxy for the rotation 
eff'ect, even in systems with blending occuring during eclipses, 
and derive values for and vsini. However, for systems with 
low eccentricity, blending of the spectral lines during eclipses is 
stronger than in VI 143 Cyg. Therefore, every systematic error 
in the tomography or in the subtraction of the foreground spec- 
trum due to the parameters used in the subtraction, will have a 
substantial eff'ect on the value of /3 derived. 

The second method presented in this paper, the modeling of 
the shape of the broadening functions, avoids the problems of the 
first method by taking the influence of the eclipsing star into ac- 
count explicitly. This makes the method more suitable for eclips- 
ing binary systems with low eccentricity (higher blending of the 
spectral lines). The use of the complete BFs, instead of the center 
of gravity or the Gaussian fit to the line, has the additional advan- 
tage that it makes the fitting of the parameters and vsini more 
precise. However, the derived values for these parameters are not 
necessarily free of systematic errors. Further work is needed here 
to make it possible to derive information about velocity fields on 
the stellar surface represented by macro-turbulence, or diff'eren- 
tial rotation. If diff'erential rotation is present it might be possible 
to derive some information about the inclination of the rotation 
axes towards the observer. 

With both methods we determined orbital and stellar param- 
eters which agree within their uncertainties with values derived 
in earlier studies, with the exception of the semi-amplitude of 
the secondary. The values derived for the angle agree with 
each other. We found that the stellar rotation axes in VI 143 Cyg 
are normal to the orbital plane for both components, which 
leaves the theoretical predictions for the apsidal motion un- 
changed. This makes V1143Cyg the first binary system with 
main-sequence stars for which the orientations of the rotation 
axes of both components are determined. 

Even with progress made in the field of apsidal motion in the 
last years (e.g. Claret & Willems 2002), it is interesting to apply 
the methods discussed here to other binary systems where there 
is a stronger disagreement between the observed and expected 
apsidal motion (e.g. DI Herculis, Claret 1998), in order to in- 
vestigate whether a misalignment of the stellar rotation axes can 
be excluded as a cause for the disagreement, or a misalignment 
contributes to the diff'erence. Furthermore, information about the 
orientation of the stellar rotation axes in binary systems could 
be used to explore the dependence of the synchronization time 
scales on the semi-major axis, eccentricity of the binary system 
and the stellar type of its components. This might lead to new 
insight about the formation and evolution of binary systems and 
stellar interiors. 
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Table A.l. Radial velocity measurements of VI 143 Cyg out of eclipse and during the primary and secondary eclipse. The photon midpoints of the 
observations are given in the first column, and the radial velocities of the primary and secondary are given in the second and third column. 



HJD 


Vel.P 


Vel. S 


[JD-2400000] 


[km/s] 


[km/s] 


53612.665 


35.40+0.13 


-81.71+0.22 


53612.683 


30.58+0.14 


-76.43+0.23 


53651.654 


-71.40+0.14 


27.82+0.22 


53613.663 


-76.90+0.14 


33.40+0.23 


53613.771 


-78.44+0.14 


34.56+0.22 


53613.880 


-79.27+0.14 


35.48+0.23 


53934.843 


-79.49+0.14 


35.83+0.23 


53583.676 


-78.86+0.14 


35.39+0.23 


53936.818 


-57.40+0.13 


13.47+0.22 


53585.669 


-51.69+0.15 


6.91+0.25 


53649.639 


43.23+0.14 


-89.43+0.23 


53611.709 


62.33+0.14 


-109.40+0.23 


53611.728 


63.65+0.14 


-111.21+0.24 


53611.814 


70.69+0.13 


-117.65+0.22 


53611.835 


72.30+0.14 


-119.47+0.23 


53611.884 


76.29+0.14 


-123.12+0.23 


53611.903 


77.69+0.13 


-125.03+0.22 


53932.846 


80.45+0.13 


-127.61+0.22 


53581.675 


96.71+0.14 


-144.12+0.24 


53650.612 


88.17+0.14 


-135.32+0.23 


53650.737 


67.50+0.14 


-113.86+0.22 


53612.865 


-12.10+0.18 


-33.23+0.27 


53612.884 


-15.74+0.19 


-29.58+0.28 


53612.901 


-19.29+0.18 


-26.26+0.28 


53612.919 


-22.44+0.18 


-22.26+0.27 


53579.678 


-11.09+0.18 


-33.60+0.28 


53612.701 


27.01+0.20 




53612.719 


24.88+0.20 




53612.737 


22.93+0.20 




53612.757 


18.28+0.20 




53612.775 


6.63+0.20 




53612.793 


-2.84+0.20 




53612.811 


-5.44+0.20 




53612.829 


-6.75+0.20 




53612.847 


-8.71+0.20 




53587.717 




-46.88+0.32 


53587.732 




-46.88+0.32 


53587.748 




-46.85+0.32 


53587.767 




-47.20+0.32 


53587.782 




-47.74+0.32 


53587.798 




-48.19+0.32 


53587.816 




-50.23+0.32 


53587.832 




-52.14+0.32 


53587.848 




-53.93+0.32 


53587.866 




-56.31+0.32 


53587.882 




-57.96+0.32 



